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ANALYSIS OF DIFFERENTIAL ABSORPTION 
LIDAR TECHNIQUE FOR MEASUREMENTS OF ANHYDROUS 
HYDROGEN CHLORIDE FROM SOLID ROCKET MOTORS 
USING A DEUTERIUM FLUORIDE LASER 

Clayton H. Bair and Frank Allario 
Langley Research Center 

SUMMARY 

The present study evaluates an active optical technique (^fferential absorption 
lidar (DIAL)) for detecting, ranging, and quantifying the concentration of anhydrous HCl 
contained in the ground cloud emitted by solid rocket motors (SRM). Results are pre- 
sented of an experiment in which absorption coefficients of HCl were measured for 
several deuterium fluoride (DF) laser transitions demonstrating for the first time that a 
close overlap exists between the 2-1 P(3) vibrational transition of the DF laser and the 
1-0 P( 6 ) absorption line of HCl, with an absorption coefficient of 5.64 (atm-cm)“^. These 
measurements show that the DF laser can be an appropriate radiation source for detecting 
HCl in a DIAL technique. Development of a mathematical computer model to predict the 
sensitivity of DIAL for detecting anhydrous HCl in the ground cloud is outlined, and 
results that assume a commercially available DF laser as the radiation source are pre- 
sented. Elements of the computer model consist of a general development of the DIAL 
equations and calculations of the extinction and scattering coefficients used in the DIAL 
equations to estimate backscattered power. A multilayer diffusion model was used to 
predict cloud growth after stabilization in the atmosphere to obtain typical concentrations 
of molecular species in the cloud as a function of downwind distance from the launch pad 
for obtaining ranging information of the DIAL technique. 

Results of a computer model which show the received laser power and signal-to- 
noise ratio as a function of range for both column content and range-resolved measure- 
ments for different particle size distribution for the AI 2 O 3 particulates and H 2 O aerosols 
in the ground cloud are presented. It was determined that column content measurements 
could be made for cloud distances 8 km from the launch pad with center-line concentra- 
tions of approximately 2 parts per million in a Gaussian cloud with a 1.5 km width by using 
a commercially available DF laser and conventional lidar instrumentation. Similarly, it 
was determined that range-resolved measurements could be performed with a lidar 
located 5 km from a Gaussian cloud which has stabilized in the manner described by the 
multilayered diffusion model. 



INTRODUCTION 


Background 

The firing of large solid rocket motors for space shuttle boosters results in the 
generation of a ground cloud containing a combination of rocket combustion products mixed 
with atmospheric constituents. The composition of the ground cloud has been discussed in 
references 1 and 2, and among its constituents are included CO2, H2, H2O, N2, CO, HCl, 
and AI2O3. Because of the toxic and corrosive properties of anhydrous HCl and its acidic 
forms, it is desirable to understand the mechanisms of generation and depletion of anhy- 
drous HCl in the cloud. This is especially important for understanding the dynamics of 
the chemistry of acid formation in the cloud if absorption of HCl by H2O aerosols and 
AI2O3 particles is assumed. Also, for satisfying the environmental impact statement for 
the space shuttle program (July 1972) (ref. 3), the concentration of HCl in the ground cloud 
is an important parameter to be measured. Multilayer diffusion models that assess the 
potential ground-level toxic fuel hazard from rocket launches have been developed (ref. 4). 
These models apply meteorological data as input information and estimate ground level 
and cloud concentrations of the various constituents including HCl as a function of dis- 
tance from the launch pad. Although these concentrations vary with the size of the rocket 
motor and the meteorological conditions, ground level and cloud concentrations of anhy- 
drous HCl below 10 ppm are predicted after the cloud has stabilized. 

The behavior of the ground cloud after the combustion products of the rocket have 
been released is discussed in reference 2. Cloud formation occurs rapidly by turbulent 
mixing of the combustion products with the atmosphere. Subsequently, the cloud rises to 
altitudes determined by atmospheric stability and rocket motor type (approximately 2750 
to 3660 m for Titan inC), stabilizing at this altitude and drifting in a direction dictated 
by meteorological conditions. The actual size, shape, and range of the cloud are not sim- 
ply characterized, but for a Titan mC rocket motor, a typical cloud diameter of 1 km is 
observed at a range 1 to 5 km from the launch pad. Ground-based in situ samplers suffer 
from the disadvantages that knowledge of the cloud drift direction must be known prior to 
launch for accurate placement of the samplers and frequent cloud passage over water 
requires seacraft whose operation may be limited by rough sea warnings, (See ref, 5.) 
From this information, it appears that a ground-based remote- sensing technique with 
range of 1 to 5 km is needed, capable of detecting concentrations of anhydrous HCl of 0.1 
to 100 ppm with a spatial resolution of at least 150 m to obtain a three-dimensional map 
of the HCl concentration. (See ref. 6.) 
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Remote -Sensing Techniques 


Passive optical remote measurements of anhydrous HCl in the ground cloud have 
been reviewed in references 5 and 6. Passive instruments perform measurements 
either in emission or absorption and require independent radiation sources such as the 
Sun or the thermally emitting hot gas in the cloud. In the emission mode, measurements 
are limited to that period of time when the cloud temperature is above ambient. In the 
absorption mode, solar radiation must be utilized; this condition requires either that the 
trajectory of the cloud be known prior to the measurement or that solar radiation scat- 
tered from cloud constituents be monitored. In either case, quantitative information of 
the HCl concentration in the cloud requires solution of complex radiative transfer equa- 
tions. These limitations of passive techniques suggest consideration of active remote 
sensii^ techniques. 

In this study, several active techniques were initially considered and included 

(1) Raman scattering (RS) 

(2) Resonance Raman scattering (RRS) 

(3) Laser induced fluorescence (LIF) 

(4) Differential absorption lidar (DIAL) 

For the lower limits of concentration considered for our study (0.1 ppm at 5 km), 
it was concluded that only DIAL could provide the required sensitivity, specificity, and 
range required for the measurement. A comparative study of the active techniques has 
been performed in references 7 and 8. RS suffers from an extremely small scattering 
coefficient; thus, this technique is limited to ranges less than 1 km. RRS which should 
increase the magnitude of the scattering cross section has not been demonstrated for HCl 
and requires operation of lasers in the ultraviolet (UV) where atmospheric scattering and 
extinction are severe. LIF of HCl in the UV or infrared (IR) has not been demonstrated 
and although a potentially larger cross section than RS exists, the problem of quenching 
by collisions with foreign gases restricts the quantitative information that can be derived 
from a fluorescence measurement. 


DIAL Technique 

Differential absorption lidar has been suggested and demonstrated by several 
authors as an active ranging technique for detecting ambient concentrations of molecules 
such as H 2 O, NO 2 , CO, and SO 2 by using visible and UV lasers. (See refs. 9 to 15.) One 
version of the DIAL technique applied to the detection of HCl in a £olid rocket motor 
(SRM) ground cloud is illustrated in figure 1. Two laser pulses are transmitted to the 
pollutant cloud, each pulse sequentially displaced in time from the other by time inter- 
val At . Some of the laser photons penetrate the cloud and are backscattered by 
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AI2O2 particles and H2O aerosols in the cloud to the laser receiver. One pulse has a 
wavelength overlapping one of the vibrational- rotational absorption lines of the 
3 . 5 -jum band of HCl. The second pulse has a wavelength which falls in between the 

absorption lines of HCl but sufficiently close to to account for signal changes in 

the backscattered signal due to changes in backs catte ring. Photons at will not be 

absorbed by HCl, but serve as a reference signal for the backscattering coefficient. For 
this application of DIAL, a reference signal X°^^ is extremely important since the vol- 
ume backscattering is primarily produced by the distribution of large particles in the 
cloud (Mie scattering). The Mie scattering coefficient and distribution of particles cur- 
rently are not well known empirically. Therefore, any changes in backscattering due to 
AI2O2 particles and H2O aerosols in the cloud must be carefully monitored by a reference 
signal. 

In this paper, results of the study to evaluate the application of DIAL to detecting 
anhydrous HCl in SRM ground clouds are presented. Absorption coefficients of HCl 
obtained by using a deuterium fluoride (DF) laser as the radiation source are summa- 
rized. Demonstration that the DF laser is a suitable radiation source for this problem 
is important since the DF laser has demonstrated high energy (>1 joule/pulse) and is 
commercially available for immediate application to a field experiment with slight modi- 
fication. The development of a mathematical computer model to predict the applicability 
of DIAL to detection of anhydrous HCl is presented. The concentration and spatial distri- 
bution of ground cloud constituents as a function of time are described by a multilayer 
diffusion model. Empirical data are used to obtain analytically the apparent index of 
refraction of the cloud due to the H2O aerosols and AI2O2 particles emitted by the rocket 
motor. The variation of laser signal return and S/N ratio for a variety of cloud condi- 
tions is presented, including specific limitations imposed on the DIAL technique by large 
concentrations of H2O aerosols in the cloud. The analytical studies are summarized and 
recommendations for implementing this technique at on-site rocket launches are made. 

SYMBOLS 

A lidar receiver telescope area 

^n’^n complex scattering coefficients 

c speed of light 

Af detector bandwidth 

ij^(x,m,e) first Mie angular dependent intensity function 
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i2(x,m,e) second Mie angular dependent intensity function 
Kabs(^>^“) Mie absorption efficiency factor 
extinction efficiency factor 
K_„„(r,X,m) Mie scattering efficiency factor 
k absorption coefficient 

k(R) extinction coefficient due to all aerosol and molecular constituents other 

than pollutant 


kcaic calculated absorption coefficient 


k 


exp 


experimental absorption coefficient 




aerosol extinction coefficient 


k»(R) pollutant absorption 

m aerosol index of refraction 


NEP noise equivalent power 

Nq particle number density 

Np(R) pollutant concentration at a distance R from lidar 
average pollutant concentration in ith range cell 
AN difference between total burden measurements made at successive range cells 

n(r) particle size distribution function 

Pj any returned power term in retrieval concentration expressions 

APj. uncertainty (noise) associated with Pj 
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output laser power 


backscattered power returned to lidar receiver 

backscattered power returned to lidar receiver at reference wavelength 

p°” backscattered power returned to lidar receiver at signal wavelength 

P°”(i),P°^^(i) backscattered power returned to lidar receiver at reference and signal 
wavelengths, respectively, from ith range cell 

Pj total pressure 

Re( ) real part 

distance from lidar to ith range cell 

AR size of range cell 

r particle radius 

Sj(x,m,e) first Mie complex amplitude function 
S 2 (x,m,e) second Mie complex amplitude function 
S/N signal-to-noise ratio 

X retrieved concentration interpreted as signal in signal-to-noise expression 

AX noise associated with retrieved concentration 

X Mie size parameter 

j3(R) total backscatter coefficient from a volume element a distance R from lidar 

iS(0 = 180°),/3jgQ volume backscatter coefficient 
rj optical efficiency 
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I 


X DF laser wavelength 

V HCl absorption line frequency 

laser frequencies 

functions related to Legendre polynomials and their derivatives used in Mie 
calculations 

p ■^^2^3 ^2^ mass density 

Pq aerosol mass concentration 

aerosol absorption cross section 
aerosol extinction cross section 
pollutant absorption cross section 

ir 

asca(^A»>^) aerosol scattering cross section 
r laser pulse length 

functions related to Legendre polynomials and their derivatives used in Mie 
calculations 

At time interval between successive laser pulses 

Superscripts: 

on refers to frequency overlapping absorption line 

off refers to frequency not overlapping absorption line 

Subscripts: 

max maximum 

min minimum 
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MEASUREMENTS OF HCl ABSORPTION COEFFICIENTS 


A primary requirement for implementation of the. DIAL technique for atmospheric 
pollution sensing is the availability of a high energy pulsed laser source with a wavelength 
overlapping a vibrational -rotational transition of the molecule to be detected. HCl has a 
fundamental absorption band at 3.5 jum, and potential laser sources for DIAL are the 
optical parametric oscillator (OPO), the HCl laser, and the DF laser. For early imple- 
mentation of the DIAL technique, only the DF laser has demonstrated sufficient energy 
per pulse for ranging capabilities beyond 1 km. 

By comparing published DF laser transitions (refs. 16 and 17) with published 
values of HCl absorption line parameters (ref. 18), it appeared that several strong over- 
laps existed within the spectral region of the HCl fundamental. In order to verify the 
degree of overlap and to obtain accurate values of the HCl absorption coefficients for 
these DF laser transitions, an experiment was performed which demonstrated for the 
first time that a strong overlap exists between the 2-1 P(3) vibrational transition of DF 
and the 1-0 P(6) absorption line of HCl^"^. The measured absorption coefficient was 
5.64 (atm-cm)“^. By measuring the absorption coefficient as a function of pressure 

(ref. 19), the 2-1 P(3) transition of DF was estimated to lie within 1.32 GHz of the center 
37 

of the HCl transition. The measured absorption coefficient was used in subsequent 
DIAL calculations to estimate performance of this technique for detecting HCl in the 
exhaust of the SRM. 

A schematic of the instrumentation used to measure absorption coefficients of HCl 
with a DF laser is shown in figure 2. A commercial line tunable DF laser operating in 
the fundamental transverse mode was used. Several characteristics of the laser are 
given subsequently. When operating with DF as the active medium, 24 laser lines were 
observed that ranged in vibrational transitions from 4-3 to 1-0 and covered the spec- 
trum from 3.493 to 4.100 |utm. Most DF lines exhibited energies £50 mJ/pulse with a 
measured optical pulse width at the half intensity point of 500 ns. Maximum pulse ener- 
gies were obtained at a repetition rate ^.5 Hz. The laser emission of the DF laser was 
observed to be superradiant even with the use of a line selection grating. Single wave- 
length operation was obtained by using a 0.3-m spectrograph as an optical filter. As 
indicated in figure 2, the laser output was focused by a lens system onto the entrance slit 
of the spectrograph. The output of the spectrograph was collimated with a CaF 2 lens, 
and then divided into a sapphire beam splitter. Reference detector B was used to mon- 
itor amplitude fluctuations of the laser; signal detector A was used to monitor the inten- 
sity of the laser beam transmitting through the gas cell. The absorption cell was made 
of Pyrex with sapphire windows. The cell length was nominally 1 m. Premix concen- 
trations of HCl and N 2 at a total pressure of 1 atm were used. Two mixtures with 0.27 
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and 10.2 percent HCl were used. These concentrations were dictated by the availability 
of the mixtures. 

The reference and signal detectors, A and B, respectively, were pyroelectric 
detectors with focusing lenses. To process the optical signals at a 0.5-Hz repetition 
frequency and to suppress the effect of amplitude fluctuations in the laser, a two-channel 
gated integrator detection system was employed. A schematic of the instrumentation is 
shown in figure 2. Briefly, the output of each detector was sampled by a gated integrator 
and fed to an A/B ratio unit. The output of the ratio unit could be displayed directly as 
an analog signal or converted to a digital signal with an analog-to-digital converter (ADC) 
and recorded by a line printer. Typically, the digitized technique was used where, for 
each measurement, 100 laser shots were recorded; data corresponding to obvious mis- 
firings of the laser were discarded. The mean and standard deviations for each set of 
measurements were, computed. Typically, uncertainty of the reported data is better 
than ±5 percent and is primarily due to the amplitude fluctuations exhibited by the laser. 
The linearity of the detection system was also checked over the range from 1 to 100 mJ 
by comparing the energy of the laser recorded by the pyroelectric detectors with a 
calibrated thermopile. 

Prior to performing the measurements, a theoretical computer program was used 
to identify the potential overlap of the 24 DF laser lines with HCl absorption lines for 
conditions simulating absorption through a 1-m cell for the HCI-N 2 mixtures used in 
the experiment. For these calculations, values of DF laser transitions were obtained 
from the literature (refs. 16 and 18). Spectroscopic parameters for HCl were taken 
from reference 18 with a Voigt profile used to describe line shapes. Calculations of HCl 
absorption coefficients included wing effects from lines up to 22 cm"^ from the line 
center. 

A summary of the theoretical and e3q>erimental results is shown in table 1. The 
first two columns identify the DF transitions which potentially showed the strongest over- 
lap. The third and fourth columns list the experimentally measured absorption coeffi- 
cient (kexp) a^d the calculated absorption coefficient (kgaic) (atm-cm)"^, respectively, 
for a 10.2-percent mixture of HCl in N 2 . The fifth and sixth columns list similar data 
for a 0.27-percent mixture. The difference in the calculated values for the two mixtures 
is due to self -broadening effects. In the latter case because of the low concentration of 
HCl, only the 2-1 P(3) transition showed significant absorption m the experiment as 
predicted by the computer model. Other DF lines not listed in table 1 (19 lines) were 
examined for potential overlap with the 10.2-percent mixtures, but no overlaps were 
identified. The result is again in agreement with the predictions of the computer model. 
The agreement between the calculated and experimental absorption coefficients is very 
good. The data listed in table 1 indicated that the most likely candidate for in a 

remote sensing lidar measurement is the 2-1 P(3) DF laser line, since it has the largest 
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absorption coefficient. However, if this line is too strongly absorbed, the 2-1 P(4) DF 
laser line would provide a very good alternative for In the subsequent computer 

analysis, the 2-1 P(3) transition will be used as for all the cases considered. 

Another parameter which can readUy be obtained from this experiment is the 

37 

separation between the 2-1 P(3) DF transition and the HCl absorption line. An exper- 
iment of this type and the results obtained are described in detail in reference 20. Fig- 
ure 3 shows results of the experiment in which a mixture containing 2 torr of HCl and 
758 torr of N£ was reduced in pressure in steps to obtain the dependence of absorption 
coefficient (atm-cm)“^ on total pressure. The circles represent the experimental 
points. The top solid curve is the predicted dependence for a DF laser frequency, 

Ut = 2750.09 cm”^ (ref. 16), and the bottom solid curve is the predicted dependence for 

— 1 37 

1^2 = 2750.05 cm" (ref. 17), by assuming that the line position of the*HCl absorption 

line is at y(l-0 P(6)) = 2750.13 cm"^. As shown in reference 20, which also includes 
effects due to pressure shifting of the absorption line, the relative position of the DF 
laser frequency to the absorption line is 1.32 ± 0.15 GHz. This result is in good agree- 
ment with the separation obtained by assuming the 2-1 P(3) DF transition and 1-0 P(6) 
HCl line position given, respectively, by Long, et al. (ref. 16) and Benedict, et al. 

(ref. 18). 

')/ 

COMPUTER MODEL DEVELOPMENT 

In this section, the computer code developed to evaluate the feasibility of the DIAL 
technique for measuring concentrations of HCl in the ground cloud of the solid rocket 
motors (SRM) is outlined. The DIAL equation (ref. 14) applied to this problem is dis- 
cussed for range -resolved measurements of HCl from the edge of the cloud to some 
range R^ into the cloud. Calculations of the extinction and scattering coefficients used 
in the DIAL equation to estimate molecular concentrations and backscattered power are 
outlined. Elements of a multilayer diffusion model to describe the distribution of cloud 
constituents spatially and with changing time are then discussed (ref. 4). 

DIAL Equations 

The DIAL equation applied to the problem of detecting molecular constituents in a 
localized cloud with a uniformly mixed distribution has been discussed in reference 14 
for CO and NO 2 . In this work, the scatterers were assumed to be Rayleigh and Mie 
scatterers in a natural atmosphere. For detection of anhydrous HCl in the ground cloud 
of an SRM, a- more complex calculation is required to account for the complex composi- 
tion of the cloud ^HCl, AI 2 O 2 , and H20j, its change in composition with time as the cloud 
diffuses away from the laimch pad, and the change of composition as the leading edge of 
the cloud is penetrated. 
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In general, the returned power from a range cell located a distance Rj from 
the laser transmitter and receiver is given by 


-2C‘k(RHkp(R)cm 

Pp = AR /3(Ri) e ° (1) 

where rj is the optical efficiency of the receiver telescope, A is the area of the 
telescope, is the output power of the laser at the transmitter, AR = ct/2 is the 

length of range cell, c is the speed of light, t is the length of the laser pulse, and 
/3(R) i§ the backscatter coefficient including Rayleigh and Mie scattering from all con- 
stituents within the range cell AR. The extinction coefficient k(R) includes extinction 
from molecular and aerosol constituents other than the pollutant to be detected up to 
range R^; kp(R) includes molecular absorption due to pollutant within the laser path 
up to range Rj. The designations for AR and Rj are noted in figure 1. 

In the DIAL technique, two wavelengths are chosen with a wavelength on the absorp- 
tion line X°” giving rise to an absorption coefficient k?”(R), and with a wavelength off 

-P-P ^ r^■p■P 

the absorption line giving rise to an absorption coefficient (R). By choosing 

rv-f-f 

the proper DF transition, kp ^(R) « 0, If it is assumed that kp is the only wavelength- 
dependent term in equation (1), then the ratio of backscattered power at the two wave- 
lengths is given by 



Since Op*^ Np(R) = k^”(R) where Np(R) is the concentration of the molecule (cm"^) 
and 0 ®*^ is the absorption cross section of the on wavelength, then 



where the integral represents the total concentration of absorbing molecules from the 
edge of the cloud to a range Rj in the cloud. For a uniformly mixed cloud with scat- 
terers only in the natural atmosphere, equation (3) can be used directly to obtain results 
as in reference 14. For the ground cloud problem, appropriate values of k(R) and /3(R) 
must be evaluated. The methodology for calculating realistic values for k(R) and /3(R) 
is given in the section ’’Computer Model Development.” 
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For obtaining range -re solved measurements of HCl in the cloud (for example, a 
nonuniformly mixed cloud), the difference between the successive total burden measure- 
ments is found. Range -resolved measurements are generated by incrementing the 
penetration distance of total burden measurements into the cloud by steps of AR. Then 
equation (3) can be used to obtain the column content from the edge of the cloud to R 2 
where R 2 = Rj + AR. In this case equation (3) becomes 



The difference between total burden measurements AN is given as 


r 


AN = Np(R) dR - Np(R) dR = ^Aog, 


■Rl 


P°^^(2) 


P°”(2) 


logp 


P?^(l)| 


p2"(i) 


(5) 


or 


AN 


= 1 ; 


■Rl+AR 


Rl 


Np(R) dR 


2a°" 


logp 


P°^^(2) P°*^(l) 


pon(2) poff(i) 


( 6 ) 


where P°^^(l) and P°*^(l) refer to power measurements at Rj, and P°^^(2) and 
P^*^(2) refer to power measurements at R 2 . Assur 
the range cell AR leads to the following expression 


Assuming that N (R) is constant across 


N (Ri + ^ log, 

P\ 2 / oAR^on 


2ARa°*^ 


P°^^(2) P°^(l) 


pon(2) poff(i) 


(7) 


By varying the position of range cell 2 with respect to the leading edge of the cloud, a 
spatial profile of the HCl concentration from the edge of the cloud to a penetration dis- 
tance Rl into the cloud can be obtained. 

The signal-to-noise (S/N) ratio of the backscattered signal has been calculated as 
one of the output parameters of the computer code. In the 3.5-jum region where HCl has 
its fundamental absorption band, indium antimonide (InSb) is an appropriate detector. 
Reference 14 has adequately treated the noise terms in DIAL systems in the infrared. 
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According to reference 14, an infrared detector is dark current or background limited. 
The noise term AX is given by 



where X is the signal given by either equation (3) or (7) (i = 2 for column content mea- 
surements and i = 4 for range-resolved measurements) and APj = NEP^if indepen- 
dent of Pi where NEP represents the noise equivalent power and Af is the band- 
width. With the expression for noise given, the S/N becomes 


S. 

N 


X 

AX 


loge 




NEP\/Af , 


1 

2 

1 

r 1 1 

2/ 

\^[p?^^(i) 


po^d)^ 

] 


1/2 


for column content measurements and 


(8a) 


S _ X 
N AX 


log. 


P°^^(2) Pj. 


on 


(1) 


L.Pr (2) P°“(l) 


fr - 

2 

1 

2 

i 

1 

2 

1 

2 ) 

^[poff(2)_ 


P°’^(2) 

T 

_pO«(l)_ 

+ 

_pOi^(l)_ 

) 


1/2 


(8b) 


for range -resolved measurements. 

An error analysis for the DIAL technique has been published by Schotland (ref. 21) 
in which he has included uncertainties in the measurements associated with the absorption 
coefficient, atmospheric parameters, and jitter in laser frequency. A similar analysis 
should be performed for this problem, but the uncertainty in composition of the particle 
scatterers and uncertainty of condition in the ground cloud currently make this type of 
analysis difficult to apply. In the application to detection of H2O in the atmosphere, 
Schotland found that primary sources of uncertainty were uncertainties associated with 
the laser frequency, and that noise associated with the signal limited the accuracy at 
higher altitudes. 
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Calculation of Extinction and Backs cattering Coefficients 

The absorption and scattering mechanisms considered in this work for the incident 
DF radiation are shown schematically in figure 4. Absorption and scattering of 
and by the natural atmosphere were included over a rai^e extending from the lidar 

transmitter through the cloud. Absorption and scattering of the laser radiation by AI2O2 
particles, H2O aerosols, and HCl molecules were assumed to, be confined within the 
exhaust cloud. 

Molecular absorption due to natural atmospheric constituents and HCl from 
2500 cm”^ to 2900 cm"^ was computed by using a computer code developed by Seals 
(ref. 22). Line parameter data were obtained from the AFCRL data tape compilation 
(ref. 23) for the major absorbers in this region: HgO, O3, N2O, CH^, and CO2. Similar 
data were obtained for HCl from reference 18. A Voigt profile was generally used in the 
calculations and wing effects were included up to 25 cm”^ from the center of the DF laser 
line. Continuum absorption due to N2 and H2O was included. 

Results of the calculated atmospheric absorption coefficients for several DF laser 
lines considered as potential wavelengths are shown in table 2, together with some results 
of calculations performed by Long (ref. 16) and measurements made by Spencer (ref. 24). 
The absorption coefficients due to Spencer listed in the table include a calculated contri- 
bution from the N2 and H2O continua. 

In order to calculate the contribution to the extinction and backscatter due to H2O 
aerosols and AI2O2 particles, the complex index of refraction and particle size distribu- 
tion must be known at the wavelengths of the transmitting laser radiation. Since at 3.5 ]u,m 
the laser wavelength and particle size are approximately the same, Mie theory must be 
employed instead of small particle scattering (Rayleigh) or the geometric optics approxima- 
tions. From Mie theory (refs. 25 and 26), the extinction coefficient (cm"^) for a system of 
particles of different radii characterized by a size distribution n(r) in particles/cm^-jiim 
is given by 





(9a) 


where CTgjjj.(r,X,m) is the extinction cross section of a particle of radius r. As is 
indicated, is a function of X and m, where m, the index of refraction, is also 

dependent on X. Expressions similar to equation (9a) for absorption and scattering 
coefficients of aerosols and particulates can also be written. 

In order to calculate kg^^j. from Mie theory, it is first necessary to define the 
efficiency factor Kgj^^(r,X,m) where 
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(9b) 


Kext(r,X,m) = 




TTr^ 


‘i I‘2n^l)Re(a„*b„) 


n=l 


Similar expressions for the scattering and absorption efficiency factors are given as 






nr^ 


= ^ I ^ «(KP + W^) 

n=l 


(9c) 


H ^abs(^>^”^) ^ Kgj^(r,X,m) - Kg^^(r,X,m) 


7rr‘ 


(9d) 


In equations (9a) to (9d), m is the complex index of refraction, x = 2jrr/X is the Mie 
size parameter for a particle of radius r and for incident radiation at wavelength X. 
Re( ) indicates that the real part of the sum, + bj^ of the Mie complex scattering 
coefficients is to be taken. The complex scattering coefficients a^^ and b^^ can be 
ejqjressed as Bessel functions of the first and third kind and their derivatives. Values 
for Kgg^, and Kgjjg were calculated according to a technique suggested by 

Kattawar and Plass (ref. 27). 

Another parameter that must be calculated for the DIAL equation is the volume 
backscattering coefficient /3 (R) which was calculated from the angular dependent 
intensity functions ij^ and I 2 defined by 


ij^(x,m,0) = Sj(x,m,0) Si*(x,m,0) (10) 

ig(x,m,0) = S 2 (x,m, 0 ) S 2 *(x,m,e) ( 11 ) 

where 

00 

Si(x,m,0) = ^ + bnXn) (12a) 

n=l 

00 

Sj(x,m,e) = 2 (Vn »nT„) «2b) 

n=l 

and the asterisk indicates the complex conjugate. 
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In equations (12a) and (12b), and are functions related to the Legendre 
polynomials and their derivatives. Their exact dependence is given in reference 25. The 
volume backscattering coefficient at 6 = 180° for a specified particle size distribution 
and for a specified scatterer is given by 

^{e = 180°) = -^( ^ n(r) dr (13) 

47t2\ 2 / 

where ij^ and ig were calculated at 6 = 180°. The backscattering coefficient /3(R) 
used in equation ( 1 ) was calculated by adding contributions from the three scatterers, 
AlgOg, HgO, and naturally occurring particulates in the atmosphere. 

Values of the extinction and backscattering coefficients were calculated on a 
CDC 6600 computer for the five DF laser wavelengths listed in table 1. For AlgOg, four 
different particle size distributions were calculated, designated I to IV. Particle size 
distributions labeled I, n, and III were taken from the literature (refs. 28, 29, and 30, 
respectively) as suggested particle size distributions for relatively small rocket motors. 
Distribution IV was obtained from reference 31 by applying the fitted curve of recent 
experimental measurements (ref. 32) of particle size against rocket motor size to the 
Titan inC nozzle. In table 3, pertinent data for the various size distributions used in 
the calculations are shown. The first column lists the designation and appropriate refer- 
ence; the second column lists the distribution function; the third and fourth columns list 
the minimum and maximum radii assumed for the integration points, respectively; the 
fifth column lists the ratio of the mass concentration to the particle number density 
Po/^o given by 


^o _ 47rp 



(14) 


where p, the mass density of AI 2 O 3 , was taken to be 3.7 g/cm3. The leading constants 
in the distribution functions listed in the second column have been normalized so that 
Pq = 1 mg/m3. The listing in the fifth column is a convenient way to convert from liter- 
ature values given in mass density to number densities. In table 4, the values of kext 
and i 3 j 0 Q for the four AI 2 O 3 distribution functions for the five pertinent DF laser tran- 
sitions are given and were used in the DIAL equation to evaluate backscattered power 
and S/N ratio. For the AI 2 O 3 particulates, the index of refraction was taken from 
Bauer and Carlson (ref. 29). 
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A procedure similar to that used for AI2O3 was used to obtain the contribution to 
kext /3(R) for the H2O aerosols. The index of refraction for H2O aerosols was 

taken from Centeno (ref. 33). The particle size distributions marked "Cloud" and 
"Haze M" by Deirmendjian (ref. 34) were used to describe H2O aerosol size distribu- 
tion; these two models were chosen to give results representative of two vastly differ- 
ent particulate sizes. 

For the natural aerosols in the atmosphere, the indices of refraction and the size 
distribution were taken from McClatchey (ref. 35). In table 5, k^^^ and ^(R) for 
the HgO aerosols and atmospheric aerosols are given. 

Multilayer Diffusion Model 

The computer model discussed and the ejq)erimental results obtained were used 
to evaluate the feasibility of DIAL for detecting HCl in the ground cloud of solid rocket 
motors. The parameters used to describe the ground cloud are based on results of the 
multilayer diffusion model (ref. 4) developed to predict cloud growth after stabilization 
based on initial cloud dimensions and existing meteorological conditions. In general, the 
model assumes diffusion of a Gaussian plume with provisions made for vertical mixing 
and depletion due to gravitational settling and precipitation. Figure 5 shows some pre- 
liminary results from this model which have been used to provide estimates of cloud 
growth and composition as a function of range for use in the computer code for the DIAL 
calculations. The data have been obtained by assuming a Titan IOC launch vehicle during 
a sea breeze condition. Figure 6 shows the width of the Gaussian cloud (measured from 
the two points at which the total mass concentration is 10 percent of the center-line con- 
centration). Typical cloud dimensions vary from <1.2 km at 1.0 km from point of sta- 
bilization to >20 km at 14 km after stabilization. 

Figure 5 shows the predicted decrease in the concentrations of the three cloud con- 
stituents considered important in this study, AlgOg aerosols, HgO aerosols, and anhy- 
drous HCl, as a fimction of downwind distance from the point of stabilization. The rela- 
tive mass concentrations are based on the fuel properties characteristic of the Titan HIC 
launch vehicle and are not dependent on the size distributions chosen for the aerosols. If 
the actual number density is required, equation (14) can be used with p and chosen 
appropriately to AlgOg or HgO depending upon the constituents desired. 

RESULTS 

In order to evaluate the feasibility of the DIAL technique, the computer code dis- 
cussed was used to obtain two types of information. First is the received power at the 
receiver as a function of range for the two DIAL wavelengths, and X°*^. (See 

eq. (1).) Second is the S/N for the received signal as a function of range for both 
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column content and range -re solved measurements given by equations (8a) and (8b), respec- 
tively. Center-line concentrations of HCl, AI2O3 particulates, and H2O aerosols for the 
Gaussian cloud were varied as a function of downwind distance as previously discussed. 
Four AI2O2 and two H2O particle size distributions were also considered. The assumed 
parameters of the lidar system are summarized in table 6. 

Figure 7(a) shows typical results of the returned power (juW) as a function of 
range (km) for the "on,” 2-1 P(3), and 'bff," 1-0 P(7), wavelengths for a cloud 2 km down- 
wind, The laser transmitter is located at zero and the edge of the cloud at 1.35 km. The 
center-line concentrations and cloud width were taken from figures 5 and 6, respectively, 
with particulate size distributions given by model IV for the Al20g particulates and cloud 
for the H2O aerosols. The decrease of the returned signal before the cloud is encoun- 
tered is due to the 1/^^ dependence, as well as to atmospheric extinction. The signal 
difference that can be observed immediately before cloud penetration is a result of a 
differential in atmospheric attenuation (see table 2) and will give a (negative) contribution 
to the retrieved HCl concentration. The sharp increase of returned signal at 1.35 km is 
due to the enhanced backscatter from the cloud particulates. After cloud penetration, 
the returned power increases until extinction (exponential in particulate concentration) 
dominates backscatter (linear in particulate concentration). Inside the cloud, the 
2-1 P(3) laser line is attenuated faster than the 1-0 P(7) line because of the presence 
of HCl which gives a signal difference that can be related to the HCl concentration by 
equation (4) for column content or equation (7) for range-resolved measurements. 

In figure 7(b), the S/N ratio defined by equation (8) is plotted as a function of 
penetration depth into the cloud. The squares represent the total burden measurements 
and the circles represent range -resolved measurements with a 75-m range resolution. 
These results indicate that for a Gaussian cloud, 2-km downwind range resolved measure- 
ments can be made up to 700 m into the cloud with a S/N > 10. 

Figures 8(a) and 8(b) are similar to figure 7 with the exception that haze M was 
used to describe the H2O aerosols. Since haze M heavily weights submicron particles, 
approximately 1200 times more particles had to be included than in Deirmendjian* s cloud 
model to give the same mass density as required by figure 5. In spite of the increased 
number of H2O aerosols, a comparison of figures 8 and 7 shows only a slight decrease 
in returned power and S/N as a function of range. 

The effect of small Al20g particles is shown in figures 9 and 10, which employ the 
Al20g model I size distribution and are otherwise identical to figures 7 and 8, respec- 
tively. As mentioned previously, the number density of Al20g particulates is adjusted to 
give the mass concentration required by the multilayer diffusion model. (See fig. 5.) 
Although the values of returned power and S/N are lower than the large Al20g particle 
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case ^Al 203 model iv), range-resolved information can still be achieved with S/N > 10 
up to 600 m into the cloud. 

In figures 7 to 10 the capability of making a measurement on a cloud that is located 
2 km from the launch site is examined; however, as the cloud continues to drift downwind, 
its size and composition change. Figures 11 and 12 indicate the effect that this downwind 
drift would have on the measurement capability of a stationary lidar located at the launch 
site. In figure 11 the maximum S/N for a total burden measurement is shown as a 
function of downwind distance, where the H 2 O aerosols are described by the cloud model. 
The separate curves correspond to the four AI 2 O 3 particle size distributions as indicated 
in the figure. As shown by the figure, column content measurements can be made on 
clouds up to 8 km downwind for the small Al 20 g particle case and on clouds up to 14 km 
downwind for the larger particles expected for exhaust engines in the Titan HIC class. 
Figure 12 is similar to figure 11 with the exception that haze M is used to describe the 
H 2 O aerosols. As before, the greatly reduced H 2 O particle size only slightly reduces 
the S/N. In all cases in figures 11 and 12 the penetration depth at which the S/N max- 
imizes varies from approximately 300 m near the launch site to approximately 1 km for 
downwind distances exceeding 8 km. Results similar to those shown in figures 11 and 12 
show that range-resolved measurements can be performed on clouds from up to 5 km 
downwind for AlgOg model I and up to 8 km downwind for AlgOg model IV. The maxi- 
mum S/N occurs at approximately the same penetration depth as for column content 
measurements. 


CONCLUDING REMARKS 

A computer model has been developed which has shown that using an experimentally 
measured absorption coefficient of HCl at the 2-1 P(3) DF laser wavelength and conserva- 
tive laser-receiver lidar parameters, HCl can be detected in the SRM tropospheric ground 
cloud. Several particulate size distributions were considered for both H 2 O and AI 2 O 3 
aerosols. In all the cases considered, it was determined that range -resolved and column 
content measurements could be made on clouds up to 5 to 8 km, respectively, from the 
launch pad with center-line concentrations ~2 ppm and width «1.5 km. Penetration depths 
of at least 500 m could be achieved with profiling information available every 75 m. The 
downwind range at which useful information could be obtained in this lidar simulation 
was limited not by lack of differential signal but rather by lack of total returned signal. 
Therefore, a higher energy laser and/or larger telescope receiver would increase the 
measurement capability over the system considered here. Also, an alternative to the 
stationary lidar considered in this study could be an aircraft-based lidar capable of fol- 
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lowing the ground cloud and therefore restricting the 1 /r 2 total power losses charac- 
teristic of the lidar signal returned from distance R. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
December 22, 1976 
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TABLE 1.- MEASURED kgxp AND CALCULATED kcalc 
ABSORPTION COEFFICIENTS IN atm-lcm-1 OF HCl 
FOR FIVE DF LASER LINES 


DF transition 

10.2% HCl in N2 

0.27% HCl in N 

cm"^ 

Identification 

^exp 

^calc 

k 

exp 

^calc 

2839.78 

1-0 P(3) 

0.067 

0.049 

— 

0.039 

2816.38 

1-0 P(4) 

.031 

.024 

— 

.020 

2750.09 

2-1 P(3) 

— 

5.00 

5.64 

5.46 

2727.30 

2-1 P(4) 

.22 

.176 

— 

.142 

2703.99 

2-1 P(5) 

.020 

.017 

— 

.014 


TABLE 2.- ATMOSPHERIC ABSORPTION COEFFICIENTS AT 
SEVERAL DF LASER WAVELENGTHS 



Long 

Spencer 

Present study 

Transition 


^calc’ 

Vy 

Ir 

‘‘exp’ 


'‘calc’ 


cm~^ 

km"^ 

cm”^ 

km'^ 

(a) 

cm"^ 

km"^ 

1-0 P(3) 

2839.78 

0.1032 


— 

2839.78 

0.1196 

1-0 P(4) 

2816.38 

.1044 


— 

2816.38 

.1062 

^2-1 P(3) 

2750.09 

.04079 

2750.0 

0.039 

2750.09 

.04071 

'^l-O P(7) 

2743.00 

.06019 

2743.0 

.058 

2743.00 

.05929 

2-1 P(4) 

2727.30 

.06773 

2727.4 

.050 

2727.30 

.07217 

2-1 P(5) 

2703.99 

.02846 

2703.9 

.033 

2703.99 

.02862 


^Absorption coefficients listed in this column include a calculated 
contribution for the Ng and HgO continua. 

bThe 2-1 P(3) and 1-0 P(7) transitions were used as and 

respectively, in the DIAL calculations. 



TABLE 3.- DESCRIPTION OF THE AI 2 O 3 AEROSOL MODELS 
USED IN THIS STUDY 



Distribution 

function 



Po/%> 

Designation 
and reference 

rmin> 

pim 

/im 

mg/m3 

particle/ cm3 

I (ref. 28) 

1115r2e-2.4r3 . 

0.01 

2 

0.006456 

n (ref. 29) 

342.7r0-5e-2.5r 

.01 

6 

.01302 

m (ref. 30) 

48.04re-l-'^0’^ 

.01 

10 

.06595 

IV (ref. 31) 

1.920e-0.65r 

.01 

20 

.3385 


TABLE 4.- EXTINCTION COEFFiaENTS (km-l/mg/mS) AND BACKSCATTER 
COEFFICIENTS (km-lsr" Vmg/mS) FOR THE FOUR AI 2 O 3 PARTICLE 
SIZE DISTRIBUTIONS CONSIDERED IN THIS STUDY 



AI 2 O 3 model I 

AI 2 O 3 model n 

AI 2 O 3 model m 

— 

AI 2 O 3 model IV 

cm"^ 

^ext 

^180 

^ext 

^180 

CD 

^180 

^ext 

^180 

2839.78 

0.2707 

0.003950 

0.4015 

0.02210 

0.2756 

0.03532 

0.1164 

0.03526 

2816.38 

.2656 

.003906 

.3999 

.02166 

.2758 

.03479 

.1165 

.03493 

2750.09 

.2440 

.003707 

.3918 

.01945 

.2766 

.03242 

.1172 

.03335 

2743.00 

.2426 

.003676 

.3913 

.01932 

.2767 

.03224 

.1172 

.03318 

2727.30 

.2390 

.003670 

.3899 

.01903 

.2768 

.03199 

.1173 

.03308 

2703.99 

.2341 

.003635 

.3880 

.01863 

.2770 

.03149 

.1176 

.03278 
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TABLE 5.- EXTINCTION COEFFICIENTS AND BACKSCATTER COEFFICIENTS 
FOR NATURALLY OCCURRING AEROSOLS AND 


TWO H 2 O AEROSOL MODELS 


I', 

Natural aerosols* 

Cloud^ 

Haze 

cm~^ 

^ext 

^180 

‘‘ext 

^180 

^ext 

^180 

2839.78 

7.115 X 10"® 

1.367 X lO'”^ 

0.3106 

6.225 X 10"® 

0.7711 

1.016 X 10"^ 

2816.38 

7.050 

1.355 

.3119 

6.325 

.7358 

9.709 X 10~® 

2750.09 

6.868 

1.338 

.3164 

7.100 

.6552 

8.840 

2743.00 

6.845 

1.335 

.3169 

7.130 

.6497 

8.782 

2727.30 

6.787 

1.333 

.3183 

7.138 

.6303 

8.507 

2703.99 

6.716 

1.327 

.3205 

7.154 

.5943 

7.998 


*Extinction coefficients and backscatter coefficients given in km"l/particle/cm3 
and km~^sr~y^particle/cm3, respectively. 

^Extinction coefficients and backscatter coefficients given in km" ^mg/m^ and 
km'^sr’V^ig/m^, respectively. 


TABLE 6.- LASER/RECEIVER PARAMETERS USED IN THE DIAL CALCULATIONS 


Laser: 





Laser output, mJ . . 
Pulse width, ns . . . 

Receiver: 

NEP, . . 

Detector area, cm2 
Telescope area, cm2 
Range gate, m . . . 
Optical efficiency . . 


2-1 P(3) DF laser line 
1-0 P(7) DF laser line 

100 

. . 500 


3 X 10"12 
. . . 0.01 
. , . 1000 
. . . 75 

. . . 0.1 
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SCHEMATIC OF LASER SIGNAL 



Figure 2.- Schematic of the HCl absorption experiment and the laser signal detection 
system used to measure overlap between HCl absorption lines and DF laser lines. 
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Absorption 

coefficient, 

k, atm cm 



Figure 3.- Variation of the absorption coefficient k with total 
pressure pj, for a 0.27-percent mixture of HCl in N 2 . 
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Rgure 4.- Extinction mechanisms considered important in the computer model for 

measuring HCI in rocket exhaust cloud. 






Downwind distance, km 

Figure 5 .- Center-line concentrations for AI2O3 and H2O aerosols and 
anhydrous HCl as a function of downwind distance for a normal 
Titan inC launch during a sea breeze regime at Kennedy Space 
Center. 
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Range, km 


(a) Signal returns as a function of range from lidar package. 

Figure 7,- Calculated signal returns and signal-to-noise ratios for a cloud located 
2 km downwind using AI2O3 model IV and H2O model cloud. 
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(b) Signal-to-noise ratios as a function of penetration depth. 
Figure 7,- Concluded. 
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(a) Signal returns as a function of range from lidar package. 

Figure 8.- Calculated signal returns and signal-to-noise ratios for a cloud located 
2 km downwind using AI2O3 model IV and H2O model haze M. 
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tO"iioiS6 r 3 ,tios 3.S 2- function of ponotmtion depth* 
Figure 8.- Concluded. 


36 


II 



Range, km 

(a) Signal returns as a function of range from lidar package. 

Figure 9 .- Calculated signal returns and signal-to-noise ratios for a cloud located 
2 km downwind using AI2O3 model I and H2O model cloud. 
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Penetration depth, km 

(b) Signal-to-noise ratios as a function of penetration depth. 
Figure 9.- Concluded. 
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(b) Signal-to-noise ratios as a function of penetration depth. 
Figure 10.- Concluded. 
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Downwind distance, km 


Figure 12.- Maximum S/N as a function of downwind distance for a 
sea breeze meteorological condition at Kennedy Space Center using 
H 2 O model haze M for four AI 2 O 3 models. 
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